Quasi-Long-Range Order in the Calogero-Sutherland Model 
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The occurrence of quasi-long-range positional order in the ground-state of the one-dimensional 
repulsive Calogero-Sutherland model is studied. By mapping the exact ground-state into a one 
dimensional classical system of interacting particles at finite temperatures the structure function 
and the displacement correlation functions are calculated numerically using Monte Carlo simulation 
methods. These are found to exhibit quasi-long-range positional order for all values of the param- 
eters. The exponent characterizing the algebraic decay of the displacement correlation functions 
with distance is estimated. It is argued that the ground-state of the repulsive Calogero-Sutherland 
model consists of a single normal phase with quasi-long-range positional order. 
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I. INTRODUCTION 

It is well established that in a classical two-dimensional 
(2D) system of particles, interacting through a potential 
that falls off sufficiently fast with distance, long-range po- 
sitional order cannot exist at finite temperatures In- 
stead, a low temperature phase with quasi-long-range po- 
sitional order (QLRPO) occurs. This phase is character- 
ized by particle displacement correlation functions that 
decay algebraically with distance, with a temperature- 
dependent exponent. This leads to a structure function 
with power law singularities at some reciprocal lattice 
vectors. It is predicted that this phase is destroyed, at a 
temperature T c , by a phase transition that can be either 
first order or continuous. The latter resulting from the 
unbiding of topological excitations. It is also known that 
theoretical considerations that take into account only 
thermal phonon fluctuations in the harmonic approxima- 
tion correctly predict the low temperature properties of 
this phase M . These predictions are confirmed by Monte 
Carlo simulations of model systems ||. 

For a quantum system of bosons in one-dimension 
(ID), interacting through a potential that falls off suf- 
ficiently fast at large distances, long-range positional or- 
der is also ruled out in the ground-state. For this system, 
theoretical analysis that accounts for zero-point phonon 
fluctuations in the harmonic approximation make pre- 
dictions analogous to those for classical 2D systems: de- 
struction of long-range positional order, power-law de- 
cay with distance of displacement correlation functions 
and power-law singularities in the structure function at 
some reciprocal lattice vectors. However, QLRPO in the 
ground-state of such ID bosons is not as well established 
as in classical 2D system, 

In this paper we study in detail positional order in 



a ID system of bosons interacting through a repulsive 
potential that varies inversely with the square of the dis- 
tance - the Calogero-Sutherland (CS) model - for which 
the exact ground-state wavefunction is known [Q. The 
structure function for the CS model was calculated in the 
harmonic approximation by several authors || . The re- 
sults confirm that it has power law singularities at recip- 
rocal lattice vectors with an exponent that varies contin- 
uously with the model parameters. However, the range 
of validity of this approximation it is not known. It is 
unclear whether QLRPO persists outside this range and, 
if so, what is the exponent governing the power-law de- 
cay of the displacement correlations. We address these 
questions here. 

In order to go beyond the harmonic approximation we 
resort to numerical calculations. Sutherland obtained 
several years ago the exact ground-state wavefunction for 
the CS model Q. Based on this wavefunction we show 
that the ground-state of the CS model can be mapped 
into a classical ID system of interacting particles at fi- 
nite temperatures. By applying Monte Carlo (MC) sim- 
ulation methods || to this classical system we calculate 
the structure function and investigate the occurrence of 
QLRPO for arbitrary values of the model parameters. 

We show evidence that the CS model ground-sate has 
a single phase with QLRPO and argue that this phase 
is non-superfluid. Over a range of parameter values we 
find that QLRPO leads to power law singularities in the 
structure function and estimate the exponents governing 
these singularities by finite size scalling analysis of the 
MC data. Outside this range we find that the structure 
function has no singularity, and that the displacement 
correlation functions decay algebraically with distance. 
We also estimate the exponent governing this decay. We 
find that this exponent differs from that predicted by the 
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harmonic approximation only in the region were quantum 
fluctuations are large, being smaller than the harmonic 
ones there. 

For particular values of the CS model parameters the 
structure function has been calculated exactly by Suther- 
land 0]. Our numerical results are found to be in good 
agreement with the exact ones. 

This paper is organized as follows. In Sec.[n] we review 
the exact results obtained by Sutherland and derive the 
mapping into a classical ID interacting system. In Sec. Ill 
we review the harmonic approximation for the CS model 
and show that Sutherland's exact ground-state wavefunc- 
tion reduces to the harmonic one in the semiclassical re- 



gion. In Sec. IV we explain our numerical method in de- 
tail and report the results obtained by it. In Sec.[v| we 
interpret these results and state our conclusions. 



II. GROUND-STATE WAVEFUNCTION 

The CS model describes N bosons of mass m on a line 
of length L interacting through the two-body potential 



where 
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According to this equation (A) may be calculated as the 
canonical ensemble average of of a ID classical 

system of particles interacting through the fictitious two- 
body potential V(x). The parameter plays the role of 
the inverse temperature. 

The numerical calculations carried out in this paper 



are based on Eq. (M), as discussed in Sec.tV 



V{r)= 1 J2 (r + nL)- 2 = ^r[sin^] 
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III. HARMONIC APPROXIMATION 



where 7 is a constant that we assume positive. This po- 
tential is periodic, with period L. For large L the n = 
term in Eq. (Q) is the most important, so that V(r) varies 
essentially as r~ 2 Q. 

Sutherland found that the exact ground-state wave- 
function for this system is given by 



const. Y\ >u < 
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where 



2A = /?=! + (! + 2g) 
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The dimensionless parameter g = -jp- measures the 
relative strengths of the potential and kinetic energies. 
Sutherland also found that the exact ground-state en- 
ergy is given by 
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E = N 
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where a = L/N is the mean interparticle distance. 

The ground-state average of any operator that depends 
only on the position operators Xi (i — 1, 2, ...,N), A({x}), 
is given by 



<^4) = 
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Using Eqs. @ and (|) | * | 2 can be written as 



Here we review the harmonic approximation for the ID 
boson system described in Sec{ll[ 



A. Phonon Spectrum and Wavefunction 

The interaction energy of the ID boson system under 
consideration here is 



(9) 



As usual, the harmonic approximation consists in ex- 
panding U to second order in the displacements, u n , from 
the classical equilibrium positions X n = na [Q. The re- 
sult is 

U = \ V[X n -X m - (u n - Um)] = E cl + SU , (10) 

n^ra 

where E c i is the classical ground-state energy and SU is 
the harmonic correction 

m = \ E d ~l^r L \*=*n-* m K - u m f . (11) 



In terms of Fourier transforms 



(12) 
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where u k is the Fourier transform of the displacement u n , 
—ir/a < k < 7r /a, and the phonon spectrum is given by 
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Approximating V(x) by V(x) — j/x 2 , the sums in Eq. 
( |l3| ) can be performed exactly. The result is 



uj{k) = 8 | k | (1 



k | a, 



(14) 



where s = (h / 'ma)ir a/ '<? / '2 is the sound velocity. This co- 
incides with the velocity of a compressional sound wave, 
s = ^/dp/dp, where p = —dE/dL is the pressure, 
p = m/a is the density and E is given by Eq. (||), if /3 is 
approximated as (3 ~ v^ff' which is justified for g> 1. 



The ground-state wavefunction in the harmonic ap- 
proximation is such that 
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Now we show that the exact wavefunction reduces to 
for j> 1. In this case (3 ~ 3> 1 and it is justified 
to apply the harmonic approximation to the fictitious po- 
tential V(x) in Eq. (^). Proceeding exactly as before we 
find 



with 



5U = \Y J ™V 2 {k)\u k | 2 , 
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where tt(k) is the fictitious potential phonon spectrum, 
given by Eq. (|l3|) with V replaced by V. Approximating 
V(x) by V(x) = — In | x |, the sums in Eq. ( |l6| ) can be 
performed exactly. The result is 



mfr(ifc) = - | k | (1 



fc | a x 



(18) 



Substituting Eq. (|T|) in Eq. (0) and using Eq. @ it 
follows that the exact wavefuntcion coincides with the 
harmonic approximation one in the limit g> 1. 



In the harmonic approximation, where the ground- 
state wavefunction, Eq. (jig), is gaussian, 



where 
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It follows from Eq. (fjj) that, 



(21) 
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Using Eq. @ we find that 
(| q(u n -u ) \ 2 ) h = 2r] h {q)[ln(2Trn) + C- Ci{2-Kn)\ 



(22) 
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where C = 0.577 is Euler's constant, Ci(x) is the cosine- 
integral function, defined as in Ref. M, and 
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The above result shows that the displacement corre- 
lation function, Eq. (|2(]), decays algebraically with dis- 
tance, that is 



(25) 



where F(n) = exp-^ /l (g)[ln (2ir) + C - Ci{2im)] is an 
oscillating function of n. From this result it follows, us- 
ing Eq. jlS|), that for q near a reciprocal lattice vector 
G p = (2-k/cl)p (p = integer) such that r) h (G p ) < 1, S(q), 
in the harmonic approximation, has power law singular- 
ities, namely 



S h (G p + K)= const. | K ^(G,)-! , 



(26) 



where K <C G t 



B. Structure Function and Correlation Function 



The structure function is defined as 



S(q) = e tqX " (e 1 ^"-^) , (19) 



where () denotes the ground-state average. 



IV. NUMERICAL CALCULATIONS 



In this Section we discuss the Monte Carlo method 
used to calculate ground-state averages of the type given 
by Eq. (ph and report the results obtained from it. 
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A. Monte Carlo Method 

Our numerical method consists of calculating averages 
from Eq. (j|) using the traditional Monte Carlo method 
P|. For a given L we first simulate the largest (3 value. 
In this case the initial configuration is an ordered chain 
of N particles. New configurations are generated by dis- 
placing particles, one at a time, in a sequential way along 
the chain. The Metropolis algorithm is used to accept or 
reject configurations. For smaller (3 values we use as the 
initial configuration the last one generated in the pre- 
vious run. Typical runs consist of 1000 MC steps per 
particle to equilibrate and 5000 MC steps per particle 
to calculate averages. We simulate systems with fixed 
a = 50 and with N ranging from N = 100 to N = 1000. 
The use of a single a value is justified because in the CS 
model the correlation functions depend only on a through 
x/a or qa. 

Our aim is to investigate by this method the behavior 
of the displacement correlation functions (e jGp ^ n "~ n °'), 
for n ^> 1. There are three possibilities. 

i- QLRPO: ( e lG p^- u ^) ->\ n p( G p) . 

ii- True long-range positional order: (e* g j>( u »» — "o)) _> 

const. . 



iii- Exponential decay: (e 



For L = oo, both true long-range positional order and 
QLRPO with 7](G p ) < 1 lead to divergencies in S(G P ). 
For finite L these singularities are replaced by peaks of 
finite height. In order to determine whether peaks corre- 
spond to singularities we perform the following finite size 
scalling analysis of the data. For a given (3 we compute 
S(G P ) for several values of L. Next we fit this data to 



S(G P ) 



const. I L 



\1-V(G P ) 



(27) 



This dependence on L arises as follows. If the system has 



QLRPO with n(G p ) < 1, Eq. (|19D predicts that S(G P ) 
diverges with L according to Eq. (27). If true long-range 
positional order is present Eq. ([l 9[) predicts S(G p ) ~ L 
and Eq. @ gives r)(G p ) = 0. 

This finite size scalling method cannot distinguish be- 
tween QLRPO with r)(G p ) > 1 and exponential decay. In 
both cases S(G P ) becomes L- independent, because there 
is no singularity. Thus, fitting the S{G P ) data to Eq. j27| ) 
results in rj(G p ) — 1. 

In order to distinguish between these two possibilities 



we compute, for a given L, 
behavior as a function of n. 



a iG p {u„~u )\ 



and study its 



B. Results 



We find that S(q), has peaks at reciprocal lattice vec- 
tors q = G p — 2irp/a, with p = 0, ±1, .., i:p max - The 



value of p m ax depends on /3 and increses with increasing 
f3. In Fig. |l| we show S(q) for (3 — 18, where two peaks 
are observed at q — G\ and q — G*2. As (3 decreases, the 
peak at q = G2 is no longer observed for /3 < 12. Below 
this value there is only one peak at q = G\. This peak 
disappears for a value of (3 < 4. 




FIG. 1. Structure function as a function of q for (3 — 18, 
L = 5000 and N = 100. Inset shows in detail the peak at 
q = G 2 - 
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FIG. 2. Exponent r/(Gi) as a function of (3. Solid squares: 
finite size scalling of MC data. Inset: typical plot used to ob- 
tain r](Gi). Continuous line: harmonic approximation. Open 



squares: ry(Gi) obtained from the decay of (e 
n for L = 5000 and N = 100. 



iGi(u„-u )\ 



with 
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To identify whether or not these peaks correspond to 
singularities we carry out the finite size scalling analysis 
described above. This is illustrated in Fig. ||. We find 
that the dependence of the exponent rj(Gi) on (3 is that 
shown in Fig. |^. 

We interpret the /3-dependence of rj(Gi) shown in Fig. 
H as follows. The numerical estimates for ?/(Gi) are in 
good agreement with harmonic approximation predic- 
tions for (3 > 5. For (3 = 4 Sutherland has obtained 
S(q) exactly. A logarithmic singularity occurs at q = G%, 
that is S(q) cx In | q — G\ |, rather than an algebraic one. 
We find that at @ = 4, r)(Gi) = 0.98 ( Fig. |). We 
believe that this is consistent with Sutherland's exact re- 
sult. The reason is that the logarithm divergence leads 
to S(G\) ~ InL in a finite system. It is well known that 
if this function is fitted to Eq. ( j27|) on a log-log plot it 
leads to an exponent equal to zero |9|] , which corresponds 
to rj(Gi) = 1. Our estimate is, within the accuracy of 
our simulations, consistent with that. For (3 — A our nu- 
merical results for S(q) agree well with the exact ones, 
as shown in Fig. || For 2 < /3 < 4 our results fitted to Eq. 
( p7| ) give i](Gi) ~ 1. As discussed in Sec.lVA, this only 



indicates that S(q) has no singularity at q = G\, which 
is consistent with either QLRPO with rj(Gi) > 1 or with 
exponential decay. 



tions. Our results for (3 — 8.0, 3.5 and 3, shown in 
Figs.2 and 4 require 10 6 MC steps. For f3 = 3 and 3.5 
we find that (e lGl '"""" ') decays with n slower than the 
harmonic approximation predicts and even becomes neg- 
ative at n = 1, as shown in Fig. 0. We also find that 
for (3 = 3.0 and 3.5 it decays as n 5 and n -10 , respec- 
tively. We interpret this as evidence that QLRPO with 
r](Gi) > I is also present for 2 < j3 < 4. Thus, there is no 
phase with exponential decay of positional correlations in 
the CS model. 
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FIG. 3. Structure function as a function of q. Compari- 
son between MC results for L = 5000 and N = 100 (solid 
squares) and exact results (continuous line) for /3 = 4 and 
(3 = 2 (inset). 

In order to investigate which of these possibilities oc- 
cur in the region 2 < [3 < 4 we compute the displace- 
ment correlation function, (e lGl ^ n_u °^), as a function 
of n, for a given L. We find that, as a result of strong 
quantum fluctuations, very long MC runs are necessary 
in order to obtain reasonably accurate correlation func- 



FIG. 4. Displacement correlation function as a function of 
n. Data points: MC data for L = 5000 and TV = 100. Con- 
tinuous line: harmonic approximation. 

For (3 = 2 Sutherland has also obtained S(q) exactly. 
In Fig. ||we show that our results for j3 = 2 agree reason- 
ably well with the exact ones. We attribute the rounding 
off near q = G\ to finite size effect. 

We also calculate the exponent //(G^) by the finite size 
scalling method for a few (3 values. We find that, in the 
region where S(G2) has a singularity, ^(Gy < 1 an d 
its estimated value is in close agreement with the har- 
monic approximation one, r q h (G2)- This approximation 
predicts that r 1 h {G 2 ) < I for (3 > 17 (g > 128). We did 
not attempt to estimate rjiG^) outside the region where 
S{G2) has no singularity from the displacement correla- 
tion function ( e i G 2 («„-«o)\_ 



V. DISCUSSION 

The main conclusion of our numerical study is that the 
ground-sate of the CS model for g > or 2 < (3 < oo has 
QLRPO characterized by an exponent that varies contin- 
uously with (3. Our results also reveal that, as far as posi- 
tional correlations are concerned, three distinct parame- 
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ter regions can be identified: i) A semiclassical region for 
[3 > 5 (g > 7.5) in which ?/(Gi) ~ rj h (G\). As discussed 
in Sec. [II this is expected for large g where the exact CS 
ground-state wavefunction coincides with the harmonic 
approximation one. Our results show that for g > 7.5 the 
CS model behaves semiclassically. ii) A region of moder- 
ate quantum fluctuations, for 4 < /3 < 5 (4 < g < 7.5), 
where S(q) has one algebraic singularity at q = G\ with 
exponent r](G\) < 1. In this region anharmonic quantum 
fluctuations are important, leading to ?y(Gi) < ri h (Gi). 
iii) A strong quantum fluctuations region, for 2 < [3 < 4, 
where S(q) has no singularities, but QLRPO exists with 
r](Gi) > 1. As a result of strong quantum fluctuations 
the displacement correlation function at short distances 
becomes negative, in contrast to the harmonic approx- 
imation that predicts positive value for all n. At large 
distances it decays algebraically with distance but with 
an exponent such that f?(Gi) <C i] h (G\). 

It is instructive to estimate the amplitude of quantum 
fluctuations in the boson positions using a 'cage model' 
fiof . This model considers a single particle moving in the 
potential well produced by the other bosons, assumed 
fixed in their classical positions. In this case we find 
that, in the harmonic approximation, the boson zero- 
point mean-square displacement from its classical equi- 
librium position is (u 2 )/a 2 = 0.55/ \Jg. According to the 
discussion above, the semiclassical region corresponds to 
g > 7.5, for which yj (u 2 )/a 2 < 0.45. We propose that 
this last result be adopted as a kind of 'Lindemann crite- 
rion' to estimate the range of validity of the harmonic 
approximation for ID interacting boson models in gen- 
eral. 

An important conclusion that can be drawn from these 
results is that, as far as QLRPO is concerned, the har- 
monic approximation gives a correct qualitative picture 
for the behavior of the CS model ground-state. The dif- 
ferences that we find between the exact result and this 
approximation, although large for /? < 5, are quantita- 
tive rather than qualitative. The fact that this conclusion 
is reached in an exactly soluble model suggests that for 
other ID boson models the harmonic approximation pre- 
dictions for QLRPO are correct, at least qualitatively, 
well beyond the semiclassical region, as long as a phase 
transition does not take place. 

One possibility that cannot be studied by our method 
is the existence of a superfluid phase. We believe that 
this phase can be ruled out in the CS model because 
the bosons are impenetrable. In ID this excludes the 
possibility of a superfluid phase. The reason is that if 
one constructs the path-integral representation for the 



ground-state partition function, the superfluid density is 
related to the boson world lines winding number fluc- 
tuations pi] . In ID these fluctuations require that the 
boson world lines cross each other, which is not possible 
if the bosons are impenetrable. From this and from our 
numerical results we conclude that the ground-state of 
the repulsive CS model has only one normal phase with 
QLRPO. We believe that a large class of ID models for 
impenetrable bosons interacting through a repulsive po- 
tential falling off so fast with distance to exclude true 
long-range positional order has a similar ground-state 
phase diagram. 
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